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ABSTRACT 

We present very high energy (VHE, E > 100 GeV) 7-ray observations of the 7-ray binary system 
LS I +61°303 obtained with the MAGIC stereo system between 2009 October and 2010 January. We de- 
tect a 6.3a 7-ray signal above 400 GeV in the combined data set. The integral flux above an energy of 300 
GeV is F(E > 300GeV) = (1.4 ± 0.3 stat ± 0.4 syst ) x lO^cm^s- 1 , which corresponds to about 1.3% of 
the Crab Nebula flux in the same energy range. The orbit-averaged flux of LS I +61°303 in the orbital phase 
interval 0.6-0.7, where a maximum of the TeV flux is expected, is lower by almost an order of magnitude com- 
pared to our previous measurements between 2005 September and 2008 January. This provides evidence for a 
new low flux state in LS I +61°303. We find that the change to the low flux state cannot be solely explained by 
an increase of photon-photon absorption around the compact star. 

Subject headings: binaries: general — gamma rays: general — stars: individual (LS 1 +61 303) — X-rays: 
binaries — X-rays: individual (LS 1+61 303) 
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1. INTRODUCTION 

The LS I +61°303 system consists of a Be star and 
a compact object of still uncertain nature, either a neu- 
tron star or a black hole. Its orbital period, which is 
most precisely mea sured in radio, is 26.4960 ± 0.0028 
days (Gregory 2002). Soft X -ray outbursts modulated with 
th e same period as in the radio waveband were reported 
by iParedes et al.l (119971) and changes in the o rbital evolution 
have been recently studied dTorres et alj|2010l) . Many other 
orbital parameters of the system are less pre cisely known 
and d i fferent solutions have been proposed (seelCasares et alj 
l2005t iGrundstrom et alj|2007t lAragona et all 120091) but ob- 
servations indicate a highly eccentric orbit (e = 0.55 ± 
0.05) with the periastron p assage at orbital phase pcr = 
0.275 (Aragon a et alj 120091) . These orbital parameters are 
important for mo delling the VHE emission of t h e sys tem 
as shown ine.g.. rSierpowska-Bartosik and Torres! d2009l) or 
iDubus et all d2010b . 

In 2006 the MAGIC collaboratio n discovered variab le VHE 
7-ray emission fromLS I +61°303 dAlbert et alJ20 06). A fol- 
lowing extensive observational campaign in Fall 2006 found a 
perio d for the VHE emission of 26.6 ± 0.2 days (Alb ert et alj 
2009). The VHE 7-ray emission shows an outburst in the 
orbital phase interval 0.6-0.7 with no significant 7-ray emis- 
sion detected during the rest of the orbit. In particular, no 
VHE 7-ray signal was detected by MAGIC around the peri- 
astron passage of the system. The data from Fall 2006 also 
sugge sted a correlati on between the X-ray and VHE 7-ray 
flux dAlbert et al.l2008l) . An extensive multi-wavelength cam- 
paign conducted in 2007, including MAGIC, XMM-Newton 
and Swift, provided strong evidence for the X -ray/VHE 7-ray 
flux c orrelation in strictly simultaneous data ( Anderhub et al. 
2009). In contrast, no correlation was found between the ra- 
dio wavelength flux and the VHE 7 -ray flux from the Fall of 
2006 campaign dAlbert et alj|2008l) . 

The VHE emission of LS I +61° 303 was confirmed 
by VERITAS observatio ns between 2006 Sep and 2007 
Feb dAcciari et alJ |2008). However, in observations con- 
ducted by the VERITAS collaboration in Fall 2008 and early 
2009, no VHE signal was detected. More recent VERITAS 
observations in Fall 2009 (the same time period as in the 
present paper) also yielde d only upper limits for VHE emis- 
sion from LS I +6T303 dAcciari et al.l 1201 lh . Very recently 
the VERITAS collaboration reported a d etection of the sys- 
tem with a significance of more than 5a dAcciari et alJl201 lh 
between orbital phases 0.05 and 0.23. This places the detec- 
tion at superior conjunction and 5.8 to 1.3 days before the 
periastron passage. No VHE 7-ray emission was previously 
detected in this phase range. 

The binary system was observe d in high energy (HE, 
0.1 - lOOGeV) 7-1-ays by EGRET dKniffen et al.lfT997h but 
the large position uncertainty of the source and inconclu- 
sive variability studies of the emission, prevented its un- 
ambiguous identification. The positional association with 
LS I +61°303 was only achieved follow ing HE 7-ray ob- 
servations by AGILE dPittori et al.l 120091) . More recently, 
FermifLP^L found that the HE 7-rays are periodically mod- 
ulate d in very good agreement with the (radio) orbital pe- 
riod (fAbdo et al. 2009) establishing beyond doubt that the sig- 
nal origins from LS I +61° 303. The HE outburst was not, 
however, observed at the same phases as the VHE outburst but 
occurred between phase 0.3 and 0.45 just after the periastron 
passage. This difference in phase may indicate that different 



processes are responsible for the HE and VHE 7-ray emis- 
sion. On the other hand, the same process might produce both 
emissions if the GeV 7-rays are produced by inverse Comp- 
ton (IC) pairs cascading developing in the radiation field of 
the star. This cascade would reduce the TeV emission and 
enhance the GeV emission when the compact obje ct is close 
to the s tar. For more details on such a scenario see iBednarekl 
(2006b). Another possibility is that the shift in the peak emis- 
sion could be caused by a different location of the 7-ray pro- 
duction site in the system (Za balza et al.ll20l"il) . We note that 
no simultaneous VHE observations are available at the same 
epoch (Au g 2008 to Jan 200 9) of the first reported Fermi ob- 
servations (Abd o et al.l2009t) . An unambiguous interpretation 
of the non simultaneous SED from MeV to TeV energies of 
LS I +61° 303 is not possible because the system might have 
changed its VHE emission in the meantime. 

Two principal scenarios have been proposed to explain the 
non-thermal emission fromLS I +61°303: an accretion pow- 
ered microquasar (e.g. iRomero et a l. 2005; Bednarek 2006a; 
iGupta and Bottchen l2006t iBo sch-Ram on et alj 120061) and a 
rotation-powered compact pulsar wind (e.g. IDubus] 120061 ; 
Sierpowska-Bartosik and Torres 2009; Zdziarski et al. 2010). 
An alternative model assumes that the compact object is an 
accreting magnetar and that the 7-r ays are produced along 
the accretion flow onto the magnetar (Bednarek 2009). High 
resolution radio measurements (Dhawa n et alj|2006 ) show an 
extended structure varying in shape and position as a function 
of the orbital phase. While this was taken as evidence for a 
pulsar wind interacting with tha t of the Be star other interpre- 
tations were suggested as well (Romero et al. 2007). Neither 
of the two proposed scenarios could be validated by accretion 
disk features, e.g. a thermal component in the X-ray spec- 
trum, or the presence of pulsed emission at any wavelength. 
Thus the engine behind the VHE emission remains an open 
question. 

Here we present new observations of LS I +61° 303 con- 
ducted with the MAGIC stereo system. This has twice the 
sensitivity of the previous MAGIC campaigns, and results in 
a significant detection of the binary system during a newly 
identified low flux state. We briefly discuss the observational 
technique and the data analysis procedure, present the VHE 
7-ray light curve of the source, and put the results in context 
of the previous VHE 7-ray observations of this system. 

2. OBSERVATIONS 

The observations were performed between 2009 Oct 15 and 
2010 Jan 22 using the MAGIC telescopes on the Canary is- 
land of La Palma (28.75°N, 17.86°W, 2225 m a.s.l.), from 
where LS I +61° 303 is observable at zenith distances above 
32°. The MAGIC stereo system consists of two imaging 
air Cherenkov telescopes, each with a 17 m diameter mirror. 
The observations were carried out in stereo mode, meaning 
only shower images which trigger simultaneously both tele- 
scopes are recorded. The stereoscopic observations provide 
a 5cr signal above 300 GeV from a source which exhibits 
0.8% of the Crab Nebula flux in 50 hours observation time, 
a factor of two more sensitive than our single telescope cam- 
paign on LS I +61°303 in 2007. Further details on the design 
a nd performance of th e MAGIC stereo system can be found 
in I Aleksic efall (1201 lh . 

The LS I +61° 303 data set spans four orbits of the sys- 
tem, with two observed for only one and three nights, re- 
spectively. The data taken in 2009 Oct and 2009 Nov were 
restricted to moonless nights. The data sample the orbital 
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phases 0.55 to 0.98 for 2009 Oct , and 0.58 to 1.02 for 2009 
Nov , the last night of which is in the next orbital cycle. The 
data recorded in 2010 Jan cover the phases 0.22 to 0.32 and 
were obtained during moonlight conditions (see Table[T]i. All 
data were taken at zenith angles between 32° and 48°. After 
pre-selection of good quality data a total of 48.4 hours of data 
remained for the analysis. The observation strategy aimed to 
cover consecutive nights with at least three hours of obser- 
vation in each individual night. Due to adverse observation 
conditions such as bad weather, the data set does not have 
uniform coverage during the orbital phases and some nights 
have shorter observation times than the planned three hours. 

3. DATA ANALYSIS 

The data analysis was performed with the standard MAGIC 
reconstruction software. The recorded shower images were 
calibrated, cleaned and used to calculated image parameters 
individually for each telescope. The energy of each event was 
then estimated using look up tables generated by Monte Carlo 
(MC) simulated 7-ray events. The events that simultaneously 
triggered both telescopes (the so-called stereo events) were 
then selectecQ and further parameters, e.g. the height of the 
shower maximum and the impact parameter from each tele- 
scope, were calculated. The gamma hadron classifications 
and reconstructions of the incoming direction of the primary 
shower particles were then perform ed using the Random For- 
est (RF) method ( I Albert et al.l2008l) . The RF calculates a vari- 
able called hadronness which is a measure of the probability 
that an event is of hadronic origin. Finally, the signal selection 
used cuts in the hadronness (calculated by the RF) and in the 
squared angular distance between the shower pointing direc- 
tion and the source position (0 2 ). The energy dependent cut 
values were determined by optimizing them on a sample of 
events recorded from the Crab Nebula under the same zenith 
angle range and similar epochs than the LS I +61° 303 data. 
For the energy spectrum and flux, the effective detector area 
was estimated by applying the same cuts used on the data sam- 
ple to a sample of MC simulated 7-rays. Finally, the energy 
spectrum was unfolded, accounting for the energy resoluti on 
and possible energy reconstruction bias dAlbert et al.ll2007l) . 

In this analysis we use for the estimation of the detection 
significance a set of cuts optimized to yield the highest sig- 
nificance on a sample of Crab Nebula data under similar ob- 
servation conditions as the LS I +61°303 data set. These cuts 
are then applied to a set of simulated MC 7-rays to estimate 
the energy threshold of the detection plot (E t h = 400GeV). 
For the light curve and spectrum determination softer cuts are 
used to reduce systematic effects and provide a lower energy 
threshold by sacrificing the highest significance. 

4. RESULTS 

The integral data set of 48.8 hours presented here results in 
a 6.3a detection of VHE 7-ray emission above 400 GeV from 
LS I +61°303 (see Fig.Q}. The integrated flux above 300 GeV 
is 

F(E > 300GeV) = (1.4±0.3 stat ±0.4 syst ) x l{r 12 cm _2 s -1 . 

corresponding to about 1.3% of the Crab Nebula flux in the 
same energy range. 

3 1 This step is only needed for the 2009 October data where no hardware 
stereo trigger was yet available. 
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FIG. 1 . — The squared angular distance between pointing direction of the 
shower and the source position (# 2 -plot) for the position of LS I +61° 303 
(points) and the simultaneous determined background regions (grey shaded 
histogram) for the total 2009/2010 MAGIC data set. N on is the number of 
events at the source position, N D g is the number of background events, N cx 
is the number of excess ev ents (N ex = N on — N Q g ) and the significance 
was calculated according to Li & Ma(1983). 

4. 1 . Light curve 

We derived a nightly light curve above an energy of 
300 GeV that is shown in Fig. [2] The measured fluxes and 
upper limits are quoted in Table Q] A constant flux fit to the 
light curve yields a x 2 /dof = 42.15/19 (p = 1.5 x 10~ 3 )and 
hence is unlikely. Thus, as in previous observations, the emis- 
sion is variable and reaches a maximum flux around orbital 
phase 0.62 of F{E > 300GeV) = (6.1 ± 1.4 stat ± 2.4 syst ) x 
10~ 12 cm _2 s _1 , corresponding to 5.4% of the Crab Nebula 
flux. This is a much lower peak emission than detected in our 
previous campaigns at the same orbital phases and sampled 
with very similar cadence. For a more quantitative compari- 
son of the 2009 emission level with the previous MAGIC ob- 
servations, we included in Figure |2]the lightcurve of the 2007 
data averaged in 0.1 phase bins. 

We found that the averaged emission level is dramati- 
cally lower than measured in our campaigns from 2005 to 
2007 dAlbert et alj|2006t lAlbert et alj 12009; An derhub et al] 
2009). Not only had the flux changed but a VHE 7-ray ex- 
cess was also observed at phases other than those of the peri- 
odic outburst between 0.6-0.7. The highest flux is, however, 
again detected in an outburst during the interval 0.6-0.7, and 
the measurements in the orbital cycle of 2009 Nov show the 
same burst profile as in previous observations but with a re- 
duced flux level. Whether the outburst recurrence is still a 
periodic property of the VHE emission for LS I +61°303 and 
whether it shows the same shape as in previous campaigns, 
cannot be determined with the small number of orbital cycles 
observed in this campaign. It is noteworthy that the outburst 
was not detected during the orbit observed in 2009 Oct. 

The mean flux for all phase bins is given in Table [2] The 
rather low mean values, even in the phase bin 0.6-0.7, of the 
individual night peak emission indicates that most of the emis- 
sion of the system is contributed by only few nights instead of 
a constant flux. A fit to a constant flux in the phase bin 0.6-0.7 
yields a x 2 /dof = 22.4/4 being strongly disfavored. 

4.2. Spectrum 

The emission level of LS I +61°303 was too low during 
most phases to obtain statistical significant phase dependent 
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FIG. 2.— VHE (E > 300GeV) 7-ray flux of LS I +61°303 as a func- 
tion of the orbital phase for the four observed orbital cycles (four upper pan- 
els) and averaged for the entire observation time (lowermost panel, black 
points). The starting MJD of each orbital cycle is given in the corresponding 
panel. I n the lower most pa nel we show as hollow triangles the previous pub- 
lished lAnderhub et al. 2009) averaged fluxes per phase bin. Vertical error 
bars show la statistical errors. 

spectra. The total signal, however, was sufficient to form a 
phase integrated spectrum with good enough statistic per bin 
to perform a Chi-square test, shown in Fig. [3] 
The spectrum is well described by a simple power law 
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with a x 2 /dof 



0.42/2. The spectral slope is compati- 
ble within errors with those previously reported by MAGIC 
(lAlbert et al]l2006t lAlbert et al.ll2009i: lAnderhub et aTll2009h . 
Hence, no evidence for long term spectral variability despite 
very different fluxes during these different campaigns is ob- 
served. Moreover, a exponential cutoff was also fitted al- 
lowing the power law parameters to vary in the one sigma 
range (adding linearly statistical and systematic uncer t ainties ) 
with respect to the fitted spectra in lAnderhub et alj (2009). 
That leads to a best fit cutoff at 483 GeV with a reduced 
X 2 /dof = 4.8/1, which is strongly increased with respect 
to the power law fit and hence strongly disfavors a cutoff in 
the spectra to explain the reduced flux level. 

5. DISCUSSION 

The binary system LS I +61° 303 was detected emitting 
VHE 7-rays in 2009-2010 at a level a factor 10 lower than 
previously observed in the phase interval 0.6-0.7. The previ- 
ously observed orbital modulated outburst in this interval was 



TABLE 1 

Observation time, orbital phase, integral flux (above 
300 GeV). Flux upper limit at the 95% co nfidence level are 
quoted in case flux significance is < 2(t jrolke et al.i2005tl . 
All errors are statistical only we estimate an additional 

systematic uncertainty of about 40%. the systematic 
uncertainty is only important in case of comparing between 
different experiments. 



Middle Time 


Obs. Time 


Phase 


Significance 


Flux 


Upper limit 


(MJD) 


(min) 




(pre-trial) 


10- 12 


10- 12 










(cm" 2 s- 1 ) 


(cm -2 s -1 ) 


55119.07 


138 


0.57 


0.8 


1.1±1.4 


4.0 


55121.08 


176 


0.64 


0.6 


0.7±1.2 


3.2 


55122.08 


194 


0.68 


1.0 


0.7±0.8 


2.4 


55126.06 


104 


0.83 


0.6 


0.8±1.3 


3.6 


55127.06 


137 


0.87 


1.9 


2.5±1.4 


5.5 


55128.07 


137 


0.91 


1.4 


1.8±1.4 


4.9 


55129.06 


132 


0.95 


3.8 


5.7±1.7 




55130.06 


135 


0.98 


2.5 


3.5±1.5 




55145.97 


140 


0.58 


2.0 


2.8±1.5 




55147.01 


221 


0.62 


4.9 


6.1±1.4 




55148.00 


216 


0.66 


-1.8 


-1.8±0.9 


1.0 


55149.02 


171 


0.70 


-0.5 


-0.5±1.1 


1.9 


55153.99 


123 


0.89 


0.5 


0.3±0.6 


1.8 


55154.98 


149 


0.92 


-0.2 


-0.2±1.2 


2.5 


55155.98 


115 


0.96 


-0.3 


-0.5±1.6 


2.9 


55156.95 


73 


1.00 


-0.7 


-1.2±1.7 


2.9 


55157.97 


82 


0.04 


0.9 


l.Oil.l 


3.6 


55215.90 


134 


0.22 


1.0 


1.3±1.4 


4.2 


55216.90 


161 


0.26 


0.6 


0.6±1.2 


3.2 


55217.90 


165 


0.30 


2.0 


2.0±1.1 





TABLE 2 

Average flux level above 300 GeV for each orbital 0. 1-phase 
bin. Flux upper limit at the 95% con fidence level are quoted 
in case flux significance is < 2(7. jrolke et al.i20051) . all 
errors are statistical only we estimate an additional 
systematic uncertainty of about 40%. the systematic 
uncertainty is only important in case of comparing between 
different experiments. 



Phase bin 



Flux 



(10- 



- 1 ) (10- 



Flux upper limit 



- 1 ) 



0.0-0.1 
0.1-0.2 
0.2-0.3 
0.3-0.4 
0.4-0.5 
0.5-0.6 
0.6-0.7 
0.7-0.8 
0.8-0.9 
0.9-1.0 



1.0 ± 1.1 

1.4 ±0.7 
0.0 ± 1.9 

1.9 ± 1.0 
1.3 ±0.5 
-0.1 ± 1.4 
1.2 ±0.7 
1.7 ±0.6 



1.3 

0.0 

1.3 

2.0 
1.1 



not detected during the first observed orbital cycle, whereas it 
was observed in the second. From this data set alone it was 
not possible to determine whether the outburst is still a truly 
periodic feature of the light curve. There are also other or- 
bital phases, varying from one orbit to another, during which 
significant emission was observed from 2005 to 2008. These 
individual nights were rarely observed and contributed only 
a minor fraction to the integral signal from LS I +61° 303 in 
these campaigns. Although we could not significantly detect 
emission at individual orbital phases in the here presented ob- 
servations because of the weakness of the source, it appears 
that several phase intervals dominate in the integral signal. 
This is strong evidence for a new behavior in the VHE 7- 
ray emission of LS I +61° 303. In previous observations con- 
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FIG. 3. — The spectrum of the complete LS I +61° 303 data sample recorded 
by the MAGIC stereo system in 2009 is shown in black. The vertical errors 
are ltr sta tistical errors. T he fit to the most precise measured spectrum from 
Fall 2006 <Albert et alJ2009l) is shown as the dashed line. The 2009 spectrum 
is fitted by a simple power law as well and the fit parameters are compatible 
to our previous measurements from 2005 to 2008 (see text for discussion). 

ducted with MAGIC, the 0.6-0.7 interval dominated the total 
flux. 

Furthermore it is evident that the flux during that phase in- 
terval is considerable reduced compared to the previous cam- 
paigns and on a similar level as in other phase intervals (e.g. 
0.9-1.0). This suggests that a change in the VHE 7-ray emis- 
sion of LS I +61°303 has occurred. On the other hand, there 
was no statistically significant change in the spectrum of the 
orbit-integrated flux in 2009 compared to the earlier results, 
suggesting that the same processes continue to produce VHE 
gamma-rays, but that either fewer are produced or they are 
more absorbed. 

If enhanced opacity causes the observed decreases in the 
VHE flux, the photons will be redistributed to lower ener- 
gies and thus might be visible as a flux enhancement. Those 
VHE 7-rays would thus need to propagate through a circum- 
source environment with a higher photon energy density and 
we would expect to detect a cut-off or an absorption feature in 
our spectrum. A cut-off in our energy range that reduces the 
integral flux by a factor of ten compare d to our previous mea- 
surements (e.g. Anderhub et al. (2009))is incompatible with 
our spectrum. The spectral fit disfavors an increased photon- 
photon absorption around the emission region as the explana- 
tion for the flux reduction in the VHE domain. 

Another possibility is that there might be fewer accelerated 
particles and/or fewer seed photons or less target matter, re- 
sulting in lower VHE 7-rays production, depending on the 



details of the assumed model (microquasar, pulsar wind, lep- 
tonic or hadronic production). However, regardless of the sce- 
nario, a change in the stellar wind density profile might ex- 
plain the change in the VHE 7-ray emission level: the wind 
density, velocity, and porosity determine the accretion rate in 
the microquasar scenario and the location of the termination 
shock in the pulsar wind scenario. Depending on the magni- 
tude of these changes it might be difficult to explain the large 
variation in the VHE domain, at least a recent study of the 
effect of wind clumping in the framework of a microquasar 
scena rio found only variations of about 10% (Owoc ki et alj 
2009). In addition, the effect of possible stellar wind density 
variations on the VHE emission in LS I +61°303 are not yet 
well understood. 

The VERITAS observations in the same period as we have 
present ed here did not dete ct VHE 7-ray emission from the 
system (lAcciari et al.11201 il) . Our measurements are not, how- 
ever, in contradiction to those of VERITAS. Our longer inte- 
gration combined with a denser sampling of two orbital cycles 
yielded a fainter detection threshold than from previous cam- 
paigns expected VHE 7-ray signal from LS I +61°303. Thus 
it is evident that a frequent sampling with long individual in- 
tegrations is required not to miss weak emission from binary 
systems. 

This is the first VHE 7-ray detection of LS I +6 1 °303 in the 
era of the Fermi satellite. The faint emission at VHE 7-rays 
does not yet permit night by night correlation studies but do 
show that the emission in LS I +6 1 °303 has changed on longer 
timescale, since 2007. More sensitive and even deeper VHE 
7-ray observations should yield shorter timescale correlation 
studies. 
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